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We study the effect of quenched coordination-number disorder of random lattices on the nature of the phase 
transition in the two-dimensional eight-state Potts model, which is of first order on regular lattices. We consider 
Poissonian random lattices of toroidal topology constructed according to the Voronoi/Delaunay prescription. 
Monte Carlo simulations yield strong evidence that the phase transition remains first order. 



1. INTRODUCTION 

Pure systems exhibiting a continuous phase 
transition are very susceptible to the addition of 
random disorder. The critical behaviour can be 
driven to new universality classes or the phase 
transition can be eliminated altogether [1]. Also 
for first-order phase transitions phenomenological 
renormalization-group arguments suggest strong 
effects caused by random disorder [2]. In partic- 
ular the order of the transition can change from 
first to second. 

The well-known paradigm to investigate such 
effects is the two-dimensional g-state Potts model 
which undergoes on regular lattices for g > 5 
a temperature driven first-order phase transition 
[3]. Monte Carlo (MC) simulations for q = 8 with 
a certain type of quenched bond randomness pro- 
vided clear evidence for a continuous phase transi- 
tion of the Ising type [4] . Also in two-dimensional 
quantum gravity studies of Potts "matter" cou- 
pled to dynamically triangulated random surfaces 
(DTRS) a similar softening effect was observed 
[5]. From a statistical mechanics viewpoint, in 
this case the Potts model is subject to annealed 
disorder in the local coordination numbers of the 
dynamical triangulation. 
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Here wc report on a study [6] of the same model 
on static Poissonian random lattices constructed 
according to the Voronoi/Delaunay prescription 
[7]. The locally varying coordination numbers 
cause the disorder similar to Ref. [5], but in our 
case the disorder is assumed to be frozen in, i.e. 
"quenched" , as in Ref. [4] . 

2. MODEL AND SIMULATION 

The 8-state Potts model is defined by the par- 
tition function 

Z=J2e-^^;E=-Y,S<r,a,;cTi = l,...,q, (1) 

(ij) 

with q = 8. The ai are integer valued spins at the 
lattice sites i, 6a-- a- denotes the usual Kronecker 
delta symbol, and the nearest-neighbor bonds (ij) 
are determined by the Voronoi/Delaunay con- 
struction of the random lattices. We always used 
periodic boimdary conditions, i.e., toroidal topol- 
ogy as depicted in Fig. 1. 

Using a standard algorithm [7,8] we generated 
20 independent replica of random lattices with 
V = 250, 500, 750, 1000, 2000, and 3000 sites and 
performed long single-cluster simulations near the 
transition point at /3 = 0.826, 0.830, 0.830, 0.830, 
0.832, and 0.833, respectively. After equilibra- 
tion we recorded 1 000 000 measurements (taken 
after 1, 1, 1, 1, 2, 4 clusters had been flipped) 
of the energy E and the magnetization M = 
{qma.x{ni} — V)/{q — l) in atime-series file, where 
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Figure 1. Random lattice with toroidal topology. 



rii < V denotes the number of spins of "orien- 
tation" i = 1,. . . ,q in one lattice configuration. 
The corresponding quantities per site will be de- 
noted by e = E/V and m = M/V. 

We then applied the reweighting method to 
compute, e.g., the specific heat, C(*'(/3) = 
((e^) — (e)^) , for each replica labeled by 
the superindex {i), performed the replica av- 

erage = = (1/20) Ef (/?), 

and finally determined the maximum, Cmax = 
C(/3cmax)- For the magnetic susceptibility, 
= I3V {[{m'^) - (m)^]), we followed exactly 
the same lines. 

The proper replica average for the specific heat 
and susceptibility follows from the general rule 
that in the quenched case the free energy (and its 
derivatives) should be averaged [9] . For the (ener- 
getic) Binder parameter, usually defined for pure 
systems as B{f3) = 1 — {e"^) /3{e^)^ , the proper 
replica average is less clear to us. We have there- 
fore studied three different definitions: Bi{fi) = 
1- [{e4)/3{e2)2], B,{I3) = 1 - [(e^)] /3 [{e^)^], 

and BaiP) = 1 - [(e")] /3 [{e^)]^ While in spin 
glass simulations [10] usually the analog of Bs 
(with e replaced by the overlap) is used, for a 
random bond Ising chain [11] a better scaling be- 
haviour was observed for the analog of Bi (with 
e replaced by m). 



3. RESULTS 

Our estimates of the extrema of C, x, and 
Bi for the various lattice sizes are collected in 
Table 1. The error bars are estimated by jack- 
kniving over the 20 replica. This takes into ac- 
count both the statistical errors on each C(')(/3) 
and the fiuctuations among the different replica. 
Already a first qualitative inspection of the data 
indicates that the first-order nature of the phase 
transition persists on quenched random lattices. 

To make this statement more precise we per- 
formed a finite-size scaling (FSS) analysis. As- 
suming a first-order phase transition, we expect 
for large system sizes an asymptotic FSS be- 
haviour of the form [12-14] 

C„,ax = ac + bcV + ..., (2) 
Xmax =ax+ + (3) 
B^,n.in = ag, + bBjV + . . . , (4) 
and 

/3c„,.. =l3t+ cc/V + (5) 

etc., where (3t is the infinite volume transition 
point. The data for C„,ax and Xm&x shown in 
Fig. 2 is clearly consistent with this assump- 
tion. From least-square fits we obtained ac = 
23.3(2.0), 6c = 0.0659(30), with a goodness-of- 
fit parameter Q = 0.16 (corresponding to a chi- 
square per degree of freedom of 1.7), and = 
-0.70(43), 6^ = 0.0629(13), with Q = 0.45. 

Also the data for the Binder parameter min- 
ima confirms the hypothesis of a first-order phase 
transition. Here the least-square fits gave asi = 
0.6240(20), be, = -18.8(1.4), Q = 0.17, as, = 
0.6236(22), fes, = -18.5(1.4), Q = 0.47, and 
as, = 0.61125(68), bs^ = -16.45(71), Q = 0.55. 
Notice the much higher accuracy of Bs . 

Our data for the pseudo-transition points and 
the corresponding fits through all data points are 
shown in Fig. 3. The resulting estimates for f3t are 
0.83360(14) from C^ax (Q = 0.51), 0.83365(14) 
from Xmax {Q = 0.47), and 0.83371(14) from 
-Bi.min (<9 = 0.40). On the scale of Fig. 3 the 
data points for I3b2„,\„ and /3B3„,i„ could hardly 
be disentangled from /3bi and are therefore 
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Table 1 

Extrema of the specific heat (Cmax), the susceptibility (Xmax), and the Binder parameter (-Bi,min), to- 



gather with the corresponding pseudo-critical couplings. 


V 




^max 




Xmax 




-^1 ,min 


250 


0.82500(44) 


33.15(45) 


0.82404(46) 


14.96(20) 


0.81872(48) 


0.5662(11) 


500 


0.82946(35) 


55.51(93) 


0.82907(34) 


31.09(56) 


0.82655(34) 


0.5875(13) 


750 


0.83087(23) 


76.1(2.0) 


0.83065(24) 


47.7(1.3) 


0.82901(24) 


0.5960(18) 


1000 


0.83112(31) 


90.4(2.6) 


0.83095(31) 


61.0(1.8) 


0.82972(32) 


0.6044(17) 


2000 


0.83232(22) 


144.8(9.0) 


0.83225(21) 


114.8(7.7) 


0.83164(21) 


0.6180(31) 


3000 


0.83300(16) 


216(11) 


0.83297(16) 


185.1(9.9) 


0.83257(16) 


0.6190(25) 



omitted. The results for /3t are 0.83350(13) from 
S2,mi„ (<3 = 0.25), and 0.83362(13) from Ba^^in 
{Q = 0.23). By taking the average of these esti- 
mates we finally obtain 

Pt = 0.83362 ± 0.00013. (6) 

Notice that this value is very close to the exactly 
known transition point of the 8-statc Potts model 
on a triangular lattice (/3j = 0.85666 . . .) [3]. 

Finally we show in Fig. 3 the "ratio-of-weights" 
definition of pseudo-transition points, Pw, which 
are expected to approach j3t exponentially fast 
with increasing lattice size [15]. Basically the idea 
is to reweight the energy histograms to a point /3w 
where the weights of the ordered and disordered 
phase are in a ratio q : 1. As in earlier studies for 
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Figure 2. FSS of specific-heat and susceptibility 
maxima. 
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Figure 3. FSS of pseudo-transition points. 

regular square lattices [14,15], we find also here 
that the /3w are quite accurate estimates of Pt 
already for very small system sizes. 

4. CONCLUSIONS 

Summarizing, we have obtained clear numeri- 
cal evidence for a first-order phase transition in 
the 8-state Potts model on quenched random lat- 
tices of Voronoi/Delaunay type. We can savely 
exclude a cross-over to a continuous transition as 
was observed for a certain type of quenched bond 
randomness on square lattices [4] and for the an- 
nealed disorder of dynamically triangulated sur- 
faces [5]. 

This conclusion is based on the FSS behaviour 
of standard thermodynamic observables. We are 
currently extending the analysis to quantities 
that are directly related to the probability distri- 
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butions of the energy or magnetization, such as 
the interface tension and the briefly mentioned 
"ratio-of-weights" definition of pseudo-transition 
points. Details of this study, which is based on a 
much larger set of 128 replica, will be presented 
elsewhere [16]. 
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